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Abstract 

We present kinetic equations that describe the evolution of 0(iV)-symmetric real 
scalar quantum fields out of thermal equilibrium in a systematic nonperturbative ap- 
proximation scheme. This description starts from the 1/iV-expansion of the 2PI ef- 
fective action to next-to-leading order, which includes scattering and memory effects. 
From this starting point one is lead to evolution equations for the propagator, which are 
nonlocal in time. Numerical solutions showed that the propagator depends only very 
smoothly on the center coordinates already after moderate times, and that correlations 
between earlier and later times are suppressed exponentially, which causes an effective 
memory loss. Exploiting these two observations, we combine a first order gradient ex- 
pansion with a Wigner transformation to derive our kinetic equations, which are local 
in time, from the nonlocal evolution equations. In contrast to standard descriptions 
based on loop expansions, our kinetic equations remain valid even for nonperturba- 
tively large fluctuations. Additionally, employing a quasi-particle approximation, we 
eventually arrive at a generalized Boltzmann equation. 



Kinetische Gleichungen aus der 
Zwei-Teilchen-Irreduziblen 1 /N-Entwicklung 

Zusammenfassung 

Wir stellen kinetische Gleichungen vor, welche die zeitliche Entwicklung von O(N)- 
symmetrischen reellen skalaren Quantenfeldern im thermischen Nichtgleichgewicht in 
einer systematischen nichtpertubativen Naherung beschreiben. Diese Beschreibung 
geht aus von der l/iV-Entwicklung der 2PI effektiven Wirkung. Die Entwicklung 
bis einschlieBlich nachstfiihrender Ordnung beriicksichtigt nichttriviale Streuprozesse 
sowie das Erinnerungsvermogen des Systems. Aus diesem Grund gelangt man zu Be- 
wegungsgleichungen fur den Propagator, welche nichtlokal in der Zeit sind. Numerische 
Ldsungen dieser Gleichungen zeigten einerseits, dafi der Propagator bereits nach mafiig 
grofien Zeiten nur noch sehr schwach von den Zentralkoordinaten abhangt, und anderer- 
seits, dafi Korrelationen zwischen friiheren und spateren Zeiten exponentiell gedampft 
werden. Diese beiden Beobachtungen erlauben es uns, durch die Kombination einer 
Gradientenentwicklung mit einer Wignertransformation zu kinetischen Gleichungen zu 
gelangen, welche lokal in der Zeit sind. Im Gegensatz zu gewohnlichen Beschreibun- 
gen, die auf einer Loop-Entwicklung beruhen, bleiben unsere Gleichungen auch dann 
noch giiltig, wenn das betrachtete System nichtpertubativ grofie Fluktuationen zeigt. 
Die zusatzliche Anwendung einer Quasiteilchennaherung fiihrt uns letztendlich zu einer 
verallgemeinerten Boltzmanngleichung. 
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Introduction and Overview 



In recent years we have witnessed an enormous increase of interest in the dy- 
namics of quantum fields out of equilibrium. Strong motivation in elementary 
particle physics comes in particular from current and upcoming relativistic heavy 
ion collision experiments at RHIC and LHC In these experiments one aims at 
the examination of the quark-gluon plasma, which is produced in a state far 
from equilibrium. In contrast to equilibrium, nonequilibrium phenomena keep 
the information about the details of the initial conditions. Therefore, their un- 
derstanding furnishes the only means to learn something about the earlier stages 
of a collision. Further nonequilibrium phenomena can be found in cosmology, 
where one is interested for example in the inflationary dynamics of the early uni- 
verse. A paradigm in this context is the phenomenon of parametric resonance in 
quantum field theory, which represents an important building block for our un- 
derstanding the (pre)heating of the early universe after a period of inflation [HI2]- 
Another very interesting problem found in astroparticle physics is the generation 
of the observed matter-antimatter asymmetry of the universe. Indeed, already in 
1967 Sakharov noticed that this asymmetry can only be generated in a universe 
out of thermal equilibrium Furthermore, closely related to the description of 
nonequilibrium dynamics there appears the fundamental question of how macro- 
scopically irreversible and dissipative behaviour can arise from microscopically 
time-reversal invariant laws of nature [HE]. Or in other words, how does a closed 
system that is initially arbitrarily far away from equilibrium equilibrate? 

There have been many different approaches to the description of nonequilib- 
rium phenomena in quantum field theory. If one can guarantee that the coupling 
is small enough, one can use perturbative descriptions based on loop expansions 
0j. This approach has its limitations in the fact that the coupling is in general 
not constant, as for example in QCD, and may become too large to be an ap- 
propriate expansion parameter. Efforts to find a nonperturbative, time-reflection 
invariant approach include the mean-field-type leading order large- N approxima- 
tion of the effective action jH]. However, this approximation neglects scattering 
and shows infinitely many conserved quantities not present in the underlying 
theory, which prohibit the approach to thermal equilibrium [5J [7j . An attempt 
to go beyond mean-field was the 1/iV-expansion of the 1PI effective action to 
next-to-leading order But this approximation can be secular (unbounded) 
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in time, and therefore fails to give a controlled description for time evolution 
problems jH]. Computations in classical statistical field theory can be used as 
an approximate means to describe the evolution of quantum fields out of equi- 
librium. These computations are expected to yield quantitatively reliable results 
if the effective particle number density is sufficiently large, such that quantum 
fluctuations are dominated by statistical fluctuations. Of course they cannot 
describe the approach to quantum thermal equilibrium |T2*1 ITS] . 

Substantial progress was achieved recently in Ref. [3], where it was shown that 
the 1 /iV-expansion of the 2PI effective action to next-to-leading order furnishes a 
controlled nonperturbative description of the dynamics of 0(iV)-symmetric real 
scalar quantum fields far from equilibrium, which is based on first principles. 
This approach does not suffer from any of the previously mentioned restrictions. 
It holds even for systems containing strongly interacting particles far from equi- 
librium. In particular, it removes spurious constants of motion encountered in 
leading order approximations and permits the description of thermalization. As 
the effective action contains all information on the quantum fluctuations of the 
system under consideration, this approach is valid even for small particle numbers 
where quantum fluctuations are dominating. These features allow one to obtain 
quantitatively reliable results even for large times. 

The starting point for this description is the 2PI effective action defined on 
the closed Schwinger-Keldysh time contour jH H2J HH1 UH UH d] • The con- 
dition, that the 2PI effective action be stationary with respect to variations of 
the propagator, leads via the Schwinger-Dyson equation to evolution equations 
for the spectral function and the symmetric propagator that are equivalent to 
the Kadanoff-Baym equations [T7J 1213 E] • These evolution equations take di- 
rect scattering into account, which leads to memory integrals. These integrals 
over time history make the evolution equations nonlocal in time. Numerical so- 
lutions of these evolution equations (see Figs. [T] and |2J) show that the evolution 
of nonequilibrium quantum fields has three generic features [Sj: 

1. effective memory loss 

(exponential suppression of unequal-time correlation functions, see Fig. |2J) 

2. intermediate-time drifting 

(smooth, parametrically slow change of modes of the equal-time correlation 
functions, see Fig. |TJ) 

3. late-time thermalization 

(exponential approach to thermal equilibrium, see Fig. 

While in equilibrium the propagator depends on the relative coordinates only, 
its additional dependence on the center coordinates is characteristic for its off- 
equilibrium behaviour. However, Fig.^shows that the evolution of the equal-time 
propagator in the intermediate-time drifting regime is very smooth. Furthermore, 
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Figure 1: Evolution of the equal-time propagator for three different 
momenta and three different initial conditions. We see intermediate- 
time drifting and late-time thermalization. Taken from jlj. 

Fig. |21 shows that already for early times correlations between earlier and later 
times are exponentially suppressed, which means that there is an effective memory 
loss. We will exploit these two properties of the propagator in order to describe 
the intermediate-time drifting and late-time thermalization by kinetic equations 
which are local in time. We obtain these kinetic equations from the nonlocal 
evolution equations by combining a first order gradient expansion with a Wigner 
transformation ^7J 1201 E] • The application of the first order gradient expansion 
is justified by the smooth dependence of the propagator on the center coordinates, 
and the effective memory loss allows us to send the limits of the memory integrals 
to ±00, which is a necessary condition for the Wigner transformation. 

The advantage of our kinetic equations is that their range of applicability is 
not limited by the validity of a loop expansion, but instead is restricted only by 
the gradient expansion that we employed in their derivation. As we will show 
explicitly, they are even valid for nonperturbatively large fluctuations |2j • In such 
a situation there are contributions with leading order in the coupling from all 2PI 
loop orders, which causes any loop expansion to break down. Compared to the 
full evolution equations, in their range of applicability our kinetic equations have 
the advantage that they do not involve an integration over time history. This is 
important for an efficient description of the late-time behaviour of nonequilibrium 
quantum fields. Apart from the description of the intermediate-time drifting 
regime, as a further very important application these equations are expected to 
be valid near a second-order phase transition. 
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Figure 2: Logarithmic plot of the evolution of the unequal-time 
spectral function and symmetric propagator. Correlations between 
earlier and later times are exponentially suppressed, which causes 
an effective memory loss. Taken from jS]. 

Additionally employing a quasi-particle approximation, we arrive at a gener- 
alized Boltzmann equation ^7J I20J ISJ I22] • Exactly as for the classical Boltzmann 
equation, its collision integral consists of terms representing the scattering of two 
quasi-particles which can be divided into gain and loss terms. However, in con- 
trast to the classical Boltzmann equation our generalized Boltzmann equation 
includes an effective mass and an effective coupling which arise from the contri- 
butions to the 1/iV-expansion of the 2PI effective action to next-to- leading order. 
Due to the effective coupling we expect our generalized Boltzmann equation to 
be valid even for large particle number densities. 

In the first chapter we review briefly the description of quantum fields in 
and out of equilibrium using the imaginary-time path and the closed Schwinger- 
Keldysh time contour, respectively [H El EH 113 CHI 121121123 In the 
second chapter we specialize to the description of far-from-equilibrium quantum 
fields using the 1/iV-expansion of the 2PI effective action to next-to-leading order. 
Here we review the derivation of the evolution equations that lead to Figs. [T] 
and El 01 IH1 12Z1 • In the third chapter we concentrate on the description of 
the intermediate-time drifting and late-time thermalization regimes and derive 
the kinetic equations and the generalized Boltzmann equation mentioned above. 
Finally, we give suggestions for a numerical implementation to solve the kinetic 
equations as well as the generalized Boltzmann equation [T7| IT^l I2TH I2T|. 



Chapter 1 
Foundations 



In this chapter we review briefly some of the rather general foundations of our 
considerations. After the introduction of our notation we consider quantum fields 
in equilibrium. The main motivation for this is Boltzmann's conjecture that any 
closed system initially arbitrarily far away from equilibrium eventually equili- 
brates, which means that its late-time behaviour can effectively be described by 
an equilibrium ensemble average. In order to check whether our kinetic equations, 
which we derive in chapter 3, contain a consistent description of the equilibrium, 
or whether the propagator has approached its equilibrium form to some given ac- 
curacy, we need to know some of the properties of the thermal propagator. These 
properties are reviewed in Sect. 1.2, using the imaginary time path introduced 
by Matsubara [23J Ell- 
in contrast to equilibrium, for a nonequilibrated system the density operator 
depends on time, and the starting point to describe such a system changes dra- 
matically. The main difference is that, while the equilibrium keeps no information 
on the details of the conditions for earlier times, the description of nonequilibrium 
phenomena is an initial value problem. However, the techniques to describe such 
a system change only slightly. One has to switch to the closed Schwinger-Keldysh 
time contour [H EG2 EE EEI1 EEE! • 



1.1 Framework and Notation 



Throughout this work we use units where Planck's constant, the speed of light 
and Boltzmann's constant are taken to be unity. We consider neutral spinless 
bosons of mass m, which can be described by a real scalar quantum field 



€ fx) 



27T 



d 3 P 



(exp (ipx) a a (p) + exp (-ipx) a\ (p)) 



The index a G {1,...,N} enumerates the different particle species, a a and a\ 
are annihilation and creation operators respectively, and p° = a/p 2 + m 2 . The 
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superscript / indicates that (x) is an interaction picture quantum field. In the 
interaction picture time translations are generated by the free Hamiltonian H : 

<f>l (x) = exp {ix°H ) $ a (x, 0) exp (-ix°H ) . (1.1.1) 

Quantum fields in the Heisenberg picture do not have any superscript at all. 
Their time dependence is governed by the full Hamiltonian H: 

$ a (x) = exp (ix°H) $ a (x, 0) exp (—ix°H) . 

The dynamics of the system under consideration are governed by the Lagrangian 
density 

Sf{x) = ± (d^ a (x) ) (>$ a (x) ) - im 2 $ a (x) $ a (x) - A ($ a (x) $ a (x) 

The plus sign of the kinetic term indicates that we use the metric where the 
time-time component is positive. Of course, we use the summation convention - 
not only for the Lorentz indices, but also for the particle species indices. Being 
a bosonic quantum field, $ a (x) satisfies the equal-time commutation relation 

[$ a (x,t);$ 6 (y,t)] = 0. 

Hence, for each t there is a simultaneous eigenstate \ip (t)) of the Heisenberg 
operators $ a (x, t) with eigenvalues cp a (x): 

$ a (x,t) \lf(t)) = V?a(x) \(p(t)) . 

The expectation value of the Heisenberg quantum field $ a (x) is denoted by 

4>a{x) = ($ (x)) = tr (V$ a (x)) . 

The trace on the right hand side represents the ensemble average with respect to 
the density operator T>, which describes the system. Accordingly, the connected 
two-point Green's function, sometimes also called connected propagator, is given 
by: 

G ab (x,y) = (T T {Qa(x)Q b (y)})-<i> a (x)<f> b (y) . (1.1.2) 

The index T indicates that the time ordering follows a certain time contour in 
the complex plane, i.e. the time coordinates are functions of the form 

z° : [0; 1] -> C ; u i-> z° (u) . 

In the course of this work we will use different time contours, the form of which 
depends on the quantity under examination. In order to perform the time or- 
dering explicitly we use functions defined on the corresponding time contour 
by 

Qt(x°M,V°(uv))=°(u*-Uv) = { I ]zZ<Z ■ 
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Similarly, we can define contour 5 functions 

d 



6 r {x° (u x ) : y° (u y )) 



(Jc0 e T (x° (u x ) y (u y )) 

-^■6 (U X - Uy) . 

dux 



However, from now on we suppress the arguments of the time variables, and think 
of them as just living on the corresponding contour being ordered according to 
their emergence on the contour. Consequently, we use the following notation for 
integrals: 



d 4 x [f (x)} 



d 3 x / du 



T 



f (x, x° (it)) 



dx° (u) 
du 



Eventually, we also generalize the definition of the functional derivative according 
to 



d 4 x 



T 



h (x) 



This definition leads immediately to 

Sf (x) 



li m -(F[f + eh]-F[f}) . 



Sf (y) 



5(x-y)5 T (x ,y )=5 T (x-y) . 



Defining the two-point correlation functions G >tab (x, y) and G < ^ ab (x, y) by 

G>,ab (x, y) = ($ a (x) (y)) - <p a (x) <p h (y) = G <M (y, x) , 

we can write the connected propagator in the form 

G ab (x, y) = r {x°, y°) G^ ab (x, y) + r {y°, x°) G^ ab (x, y) , (1.1.3) 

a decomposition that will be of great use in subsequent chapters. While the 
G functions appearing in Eq. (|1.1.3|) cause an inanalyticity of the connected 
two-point Green's function with respect to its time arguments, the correlation 
functions G >)Clb (x, y) and G <Ab (x, y) are analytic functions of their time argu- 
ments. Another very useful property of the correlation functions is their so-called 
"hermiticity" : 

G* >jab (x,y) = (tr(V$ a (x)<S> b (y))y 

= tr ($6 (y) $ a (x) V) 
= G >M (y,x) . 
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The same relation holds for G <jab (x,y). The correlation functions G >A b{x,y) 
and G <tab (x, y) are also the constituents of some other frequently used quantities, 
such as the retarded and advanced propagator 

G R , ab (x, y) = iO (x° - y°) (G >tClb (x, y) - G <Ab (x, y) ) , 

G A ,ab (x, y) = -iO (y° - x°) (G >tab (x, y) - G <>ab (x, y) ) , 
the symmetric or statistical propagator 

Fab (x, y) = ^ (G >)ab (x, y) + G <>ab (x, y) ) , 
as well as the spectral function 

Qab (x, y) = G>,a6 (x, y) - G^ab (x, y) = ( [$ a (x) ; $ fe 



We stress that we have used the usual 9 function for real time arguments in 
the definitions of the retarded and the advanced propagator. Therefore, their 
definitions are meaningful only if x° and y° are real. Because of the hermiticity 
of G >>ab (x, y) and G^ ab (x, y) the retarded and the advanced propagator are real 
functions. Furthermore, it follows immediately from the definition of G >>ab {x, y) 
and G^ah (x, y) that 

G R ,ab (x, y) = Gam (v, x ) ■ 



1.2 Thermal Quantum Field Theory 

In this section we consider a closed system in equilibrium with time-independent 
Hamiltonian H and temperature T — ^. We will work with the canonical density 
operator 1 

V = ±exp(-PH) , (1.2.1) 

where 

Z = tr(exp(-/3#)) (1.2.2) 

is the partition function. The similarity between exp {—j3H) and the time evolu- 
tion operator exp (—i (t — t ) H) suggests the definition of a new operator which 
contains these two operators as special cases. We set 

U (r) = exp (-tH) 



1 Of course, there is particle creation and particle annihilation in our system. Therefore, in 
equilibrium the particle number will be such that the free energy is minimal: 

dF(T,V,N) 
dN 

The left hand side is just the definition of the chemical potential! 
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and allow for r to take complex values. It is obvious that U (r) can be interpreted 
as an evolution operator for complex-valued times. In order to compute the 
partition function we evaluate the trace using the complete set of eigenstates of 
the quantum field $ a (x) for x° = 0: 



Z = J l[d<p a (x)(<p(Q)\U((3)\<p(Q)) . 



As indicated above, by interpreting U (r) as an evolution operator for complex 
times, we obtain the following path integral representation for the matrix element: 



(<p{o)\uvs)\<pM) 



N 



H®cp a (x) [exp(-S E [<p])] 



a=l 



</Pa(x,0) = t£a(x) 
<Pa(x,-if3) = <Pa(x) 



Here Se denotes the Euclidean action which we obtain from the classical action 
by substituting r = ix°. 



iS 



.„.]* 

o 

P 

d 3 x / dr 



1 1 A 



Thus we arrive at the following path integral representation for the partition 
function: 



N 



tpa(x,0)=tp a {x,—iP) a ~ 1 

By adding a source term to the Euclidean action, we turn the partition function 
into a generating functional: 



Z[J] 



N 



X) 



(pa(x,0) = <p a {x,—i(3) 



o=l 



/ p 

exp I -Se [(f) + J ' d 3 x J dr [J a Lp a 



Using the standard procedure one can identify the nth functional derivative of 
Z [J] with respect to J with the expectation value of a time-contour ordered 
product of n quantum fields. The time contour that we have to use here is 
the imaginary time path X shown in Fig. 11.11 For example, we obtain for the 
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Im(t) 



x o-*- 



Tn- 



4 



Re(t) 



Figure 1.1: Imaginary time path, introduced by Mat- 
subara j2H|. Although we will use different values for r , 
we will always use the same index X. 



two-point Green's function: 
1 / 5 2 Z[J] 



Z \5J a (x)8J b (y) 
1 

Z 



N 



II @<Pc (z) [fa (x) fb (y) exp (S E [<p])] 



(Tx {$ a (*) $ b (y)}) 

itr f exp (-/?#) T x {$ a (x) $ 6 (y)} 



Correspondingly, the first derivative of the generating functional with respect to 
the source gives the expectation value of the quantum field $ a (x) and we obtain 
the connected thermal propagator as in Eq. (jl.l.2j) . Alternatively, we can use the 
generating functional for connected thermal n-point Green's functions 

W[J]=]og(Z[J\) 

and obtain the connected thermal propagator as 

5 2 W [J] 



Gab (x, y) 



5J a (x)5J b (y) 



j=o 



Because of the cyclic invariance of the trace we find that for < r < (3 
(T x {$ a (x, -it) $ fc (y,0)}) 

= itr ( exp (-/3H) exp (tH) $ a (x, 0) exp (-tH) (y, 0) 

= ^tr ( exp (-PH) $ 6 (y , 0) $ Q (x, -i (r - /?)) 
= (T J {$ a (x,-z(r-/3))$ b (y,0)}) 
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and 

(f) a (x, -IT) = 4> a (X, ~l (r - (3)) . 

Hence, the connected thermal propagator satisfies the following periodicity con- 
dition: 

G ab (x, -it, y, 0) = G ab (x, -i (t-P), y, 0) . (1.2.3) 

Instead of using a complete set of eigenstates of the quantum field <3> a (x) for 
calculating tr (VTj {$ a (x) $& (y)}) we also could have used a complete set of 
eigenstates of the Hamiltonian defined by 

H \n) = E n \n) . 

With their aid we obtain 

($ a (x)$ 6 (y)) = J2(l\V\n)(n\$ a (x)\m)(m\<l> b (y)\l) 

l,m,n 

= ^i<n|$ a (x,0)|m)<m|$ 6 (y,0)|n) (1.2.4) 



X 



Z 

exp (-z (x° - y°) E m ) exp (- (/3 - i (x° - y )) E n ) 



and 



Eqs. (I1.2.4|) and (jl.2.5j) show two properties of the thermal propagator. First, we 
see that the time dependence of G >t ab {x, y) takes only x° — y° into account. The 
same also holds for G <tab (x, y) and we can conclude immediately that the ther- 
mal propagator is invariant under time translations 2 . Additionally, the thermal 
propagator also is expected to be invariant under space translations, such that we 
finally arrive at a space-time translation invariant thermal propagator. Second, 
provided the matrix elements behave well enough, the convergence of the sum in 
Eq. (jl.2.4|) is controlled by the exponential functions and we find that G >Ab (x, y) 
is an analytic function of its time arguments in the domain 

{x°-y°eC - [3 <Im(x° -y°) < o} . (1.2.6) 

The first inequality is required by the second exponential function and vice versa. 
By interchanging x and y in Eq. (|1.2.4|) G <)ab (x, y) is found to be an analytic 
function of its time arguments in the domain 

|x° - y° G C < Im (x° - y°) </?} . (1.2.7) 



2 Although in our case x° and y° are purely imaginary quantities, this statement also holds 
for real times. In that case one only has to switch to a different time contour |24l I25| . 
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In fact the domains (jl.2.6|) and (jl.2.7j) are promoted to be the domains of def- 
inition for G >t ab(x,y) and G <t ab (x,y), respectively. Knowing the domains of 
definition for G >ta b(x,y) and G <ta b(x,y) becomes important, if we look more 
closely at the periodicity condition 1)1.2.3") . We have 

G ab (x, -ir, y, 0) = G >jab (x, -ir, y, 0) , 

G ab (x, -i (r - (3) , y, 0) = G <>ab (x, -i (r - f3) , y, 0) 
and obtain the Kubo-Martin-Schwinger condition [2EJ EH] 

(x - y) = G^b (x - y, x° - y° + i(3) (1.2.8) 

by taking into account that in equilibrium the propagator is invariant under 
space-time translations. The tilde indicates the use of relative coordinates as 
arguments. When we Fourier transform the Kubo-Martin-Schwinger condition 
(|1.2.8j) with respect to the relative coordinate s = x — y, we obtain 

exp (iks) G >t ab (s) 
= exp (pk ) G <: ab (k) . (1.2.9) 
Using the Fourier transformed spectral function 

Qab (k) = G >:a b (k) - G <:a b (k) 

we can write the Fourier transformed correlation functions in the form 

G < ,ab(k) = Q ab (k)n B (k°) , (1.2.10) 

G>,ab (k) = g ab (k) (1 + n B (A; )) , (1.2.11) 
where n B (k°) is the Bose-Einstein distribution function: 

HB ^ = exp 09*6) _ I ■ 

The main results of this section are the Fourier transformed Kubo-Martin-Schwin- 
ger condition (jl.2.9J) and Eqs. ()1.2.10j) and (jl.2.11|) . which we will need in chapter 
3. We have tried to arrive at these results as quickly as possible without omitting 
important information and without giving information not necessarily needed. 
Of course, there is much more to say about thermal quantum field theory, e.g. 
one could want to calculate Green's functions with real time arguments. In this 
case one has to choose a slightly different time contour as is explained in more 
detail for example in Refs. |2*H 12*5] . 



G >t ab(k) = d A s 



1.3 Quantum Fields out of Thermal Equilibrium 
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1.3 Quantum Fields out of Thermal Equilibrium 

The striking difference between equilibrium and nonequilibrium phenomena is 
that, while the canonical density operator (jl.2.1j) is time-independent, a density 
operator describing a system out of equilibrium depends on time, 

V = V (t) 

and in general is known only for some initial time 

t ti n n = . 

Therefore, the description of nonequilibrium phenomena is an initial value prob- 
lem. One starts with a set of initial n-point Green's functions given by the initial 
density operator as 

tr(p(0)$ ai ( Xl ,0)...$ an (x n ,0)) (1.3.1) 

and follows their evolution in time. This will be our task for the rest of this work. 
In order to be able to describe general nonequilibrium initial conditions as well 
as initial thermal equilibrium, we allow for the density operator to include mixed 
states, such that 

tr (V 2 (0)) < 1 . 

In what follows we will restrict our considerations to spatially homogeneous sys- 
tems which can initially be described by a Gaussian density operator, i.e. the only 
non- vanishing initial correlations are completely determined by the one- and two- 
point Green's functions and their first-order time derivatives at the initial time. 
Hence, the set of initial n-point Green's functions we start with reads 

a (x,O) = tr(P(O)$ a (x,O)) , 
a (x,O) =tr(v(0)(d x o$ a (x))) 

\ / x°=0 

G ab (x, 0, y , 0) = tr (v (0) $ a (x, 0) $ fe (y , 0) ) - a (x, 0) fe (y , 0) , 
H ab (^0,y,0)=ti(v(0)(d x o^ a (x))^ b (y)) - <j> a (x, 0) <j> b (y, 0) , 

V / x °=y°=0 

K ab (x, 0, y, 0) = tx(v (0) (d x o$ a (x)) (d y o^ b (y)) ) - a (x, 0) 0, (y, 0) . 

V / x°=y°=0 

We want to emphasize that the restriction on Gaussian initial conditions is not an 
approximation at all. It only restricts the set of systems that can be considered. 
The treatment of higher initial correlations imposes technical complications only, 
not fundamental ones [3]. 
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Next, we face the problem of calculating an expectation value of an arbitrary 
composition of Heisenberg quantum fields for times later than the initial time. 
In order to evaluate the trace (jl.3.1j) in this more general case, we have to take 
the evolution of each Heisenberg field with respect to the initial time. Thus we 
obtain for the two-point correlation function: 

($ (x) (y)) = tr (V (0) exp (ix°H) $ a (x, 0) exp (i (y° - x°) H) 

x $ 6 (y,0)exp(-zj/°if)) 
= ti(T>(0)U (0,x°)$i(x)U (x°,y )& b (y)U (y°,0)) . (1.3.2) 

In the second line we have used Eq. (Il.l.lj) and the interaction picture evolution 
operator 

U (x°, y°) = exp (ix°H Q ) exp (-i (x° - y°) H) exp (-ty°H ) , (1.3.3) 

where we have separated the full Hamiltonian H into its free part Hq and the 
interaction V: 

H = H + V . 

Differentiating Eq. (Jl.3.3)) with respect to a; gives the differential equation 

i^- U(x°,y )=V J (x )U(x ,y°) , (1.3.4) 

where 

V 1 (x°) = exp (ix°H ) Vexp (-ix°H ) . 

The interaction picture evolution operator (jl.3.3|) is determined as the unique 
solution of the differential equation (jl.3.4|) by the initial condition 3 

U(y°,y°) = l. (1.3.5) 

The differential equation (jl.3.4j) together with the initial condition (jl.3.5|) is 
equivalent to the integral equation 

x° 

U(x°,y°)=l- i J dtiV 1 (t)U (t,y )] . 

y° 

3 This follows from the theorem of Picard-Lindelof on the existence and the uniqueness of 
solutions of ordinary differential equations. 
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By iteration of this integral equation, we obtain an expansion for U (x°, y°) in 
powers of V 1 : 



X" 

U (x°, y°) = l-ijdt [V 1 (*)] 



+ -i) 



j dt x J dt 2 [V 1 (tj) V 1 (t 2 )] +... . 



For x > y this expansion can be written in the form 



X 

U(x°,y°) = l-ij dtlV 1 



+ 



'-if 



X" X" 



Texp 



dtx / dt 2 [T {V 1 (ti)V* (t 2 )}] + ... 
\ 

/ 



-i J dt [V 1 (t)] 
\ y° 



where T denotes the usual time ordering operator along the real axis. On the 
other hand, we obtain for x° < y° 



y° y° 



+ - j dh I dt 2 



( y° 



Texp 



f {V 1 ^V 1 (t 2 )} 
\ 



+ ... 



i 

\ x° 



dt [V 1 (*)] 



/ 



Here T denotes the antitemporal ordering along the real axis, i.e. the interaction 
with the latest time argument is placed rightmost. 

Now we return to the trace (jl.3.2j) and find that all operators are sorted 
according to the order in which their time arguments appear on the time contour 
shown in Fig. 11.21 with 



max 
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t,„; =0 



max 



Figure 1.2: Closed real-time path C, introduced by Schwinger 
and Keldysh [To] . 



This is how the closed Schwinger-Keldysh time contour C comes into play in 
a very natural way when one wants to describe quantum fields out of equilib- 
rium. Consequently, we have to define the n-point Green's functions on this time 
contour. For example the connected Schwinger-Keldysh propagator is given by 



G ab (x, y) = (T c {$„ (x) $ 6 (y)}) - <p a (x) <p b (y) . (1.3.6) 

With the aid of the density operator, the generating functional for these Schwin- 
ger-Keldysh Green's functions can be written in the form 

Z v [J,K] = tr (v (0) T c { exp (i J d 4 x [J a (x) $ a {x)) 

+ \ J d ' x S dAy [$a (x) Kab (x ' y) $6 {y)] ) } 



where we have introduced the bilocal source K for reasons that will become clear 
in the next chapter, when we review the derivation of the 2PI effective action. 
Exactly as in the previous section, we evaluate the trace using the complete set 
of eigenstates of the quantum field $ a (x) for x° = 0: 



zv[j,k\ = f n^ i} (x) j n^f(y 

c,x d,y 



(0)|X?(0)|y^ (0)> 



x if 



^ (0)| T c { exp (i J d A x [J a (x) $ a Or)] 



+ l - I d 4 x I d 4 y^ a (x)K ab (x,y)^ b (y)))\\^(0)) 
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The second matrix element is given by the following path integral: 



r N r / r 

/ II ^ v?c ^ exp y \ 

. . c=l ~ 



¥> c (x, > )=</3i 1) (x) 

</3 c (x, < 0")=< ( 3i 2) (x) 



OXj) \ I I d 4 X [Sf (x) + J a (x) f a (x)] 
C 



+ % 2 j X J ^ ^ a ^ Kab ^' ^ ^ ^ ) 
c c 

The matrix element of the density operator can be written in the form 

V 1 )(0)|P(0)|^ 2 )(0))=exp( ? FM) , 
where F can be expanded in the functional sense according to 

F [<p] = a® + J otx [a« (x) <p a (x)] 
c 

a abc ( X > V> Z ) ( x ) <fb (V) <fc {z) 



+ - / d 4 x / d 4 y 



1.3.7) 



c 



c 



+ - J a*x 

c c 

+ ... . 



d y I (1 : 

c 



The values of ip a (x) at the beginning and the end of the time contour are given by 
(fa (x) and (x), respectively. Of course, the sequence of the coefficients 
contains as much information as the original density matrix. So we can rewrite 
the above trace as a functional of infinitely many nonlocal sources. However, at 
the beginning of this section we restricted our considerations to Gaussian initial 
conditions. This means for the sources that 



l ai ,...,cij 



) X jj 







(Vj > 3) . 

Even the sources and are non- vanishing only at the initial time t init = 0, 
such that we can absorb them into the external sources J and K. Anyway, we can 
only compute ratios of path integrals. In such a ratio the constant factor caused 
by cancels, and therefore, we can cancel from the functional ()1.3.7j) . 
Thus, the generating functional for Schwinger-Keldysh Green's functions is given 
by the path integral 



Zv[J,K) 



N 



[I A , 



X) 



c=l 



exp ( i j d x [if (x) + J a (x) ip a (x)} 



+ - I d 4 x j d A y [(p a (x) K ab (x, y) ip b 



'1.3.8) 
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and the expectation values needed for the calculation of the connected Schwin- 
ger-Keldysh propagator (jl.3.6J) are given by 

(Tc (*) $ b („)}> = ±j£- { f^f}, ) 

Z p [0,0J \0Ja{x)5J b {y)J J=K=Q 

and 

^ ^[0,0] V 6J a (x) ) J=K=0 • 

It is important to note that the closed Schwinger-Keldysh time contour emerges 
necessarily whenever one wants to treat nonequilibrium phenomena, regardless of 
the special circumstances (elementary particle physics, description of conductors 
and semi-conductor devices, cosmology, . . . ). 

In the next chapter we will start from the generating functional (jl.3.8|) and de- 
rive the two-particle irreducible effective action defined on the closed Schwinger- 
Keldysh time contour. After having presented the 1/iV-expansion of the 2PI 
effective action, we then are going to derive equations which determine the evo- 
lution of the Schwinger-Keldysh propagator. 



Chapter 2 



Nonlocal Dynamics of Quantum 
Fields out of Equilibrium 

In this chapter we continue our considerations on quantum fields out of ther- 
mal equilibrium. Starting from the generating functional for Schwinger-Keldysh 
Green's functions we derive the two-particle irreducible effective action defined on 
the closed Schwinger-Keldysh time contour. We review the 1/iV-expansion of the 
2PI effective action to next-to-leading order as well as the derivation of the cor- 
responding evolution equations for the propagator, which lead to Figs. [Hand m 
the introduction. As we already mentioned in the introduction, these evolution 
equations furnish a controlled nonperturbative description of the time-reversal 
invariant dynamics of quantum fields far from equilibrium. 

2.1 2PI Effective Action 

The effective action was introduced perturbatively by Jeffrey Goldstone, Abdus 
Salam and Steven Weinberg in 1962 [30]. They defined it as the sum over all one- 
particle irreducible 1 diagrams. In the following year Bryce S. DeWitt gave the 
nonperturbative definition [HI]. The physical significance of the effective action 
arises for two reasons. First, it is the generating functional for 1PI n-point Green's 
functions, and second it provides an equation of motion for the expectation value 
of the quantum field, which includes quantum corrections. Later on, this form of 
the effective action was called 1PI effective action in order to distinguish it from 
higher effective actions. 

In this section we are going to summarize briefly the basic facts about the 
2PI effective action T [(f), G}. In analogy to what we said above, it is the gen- 
erating functional for 2PI n-point Green's functions expressed in terms of the 
connected propagator G a b(x,y). This means that the diagrammatic expansion 

1 A diagram is called n-particle- irreducible (nPl), if it cannot be disconnected by cutting 
through any n internal lines. 
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of for example 



5(p a (x) 5<f) b (y) 



<i>=o 



contains only two-particle irreducible diagrams, where internal lines represent 
the connected propagator. Apart from that, the 2PI effective action provides 
equations of motion for the field expectation value a [%) and the connected 
propagator G a b (x,y), which take quantum corrections into account. 

The starting point for the derivation of an explicit expression for the 2PI 
effective action is the path integral representation of the generating functional 
for Schwinger-Keldysh Green's functions (jl.3.8|) . namely: 



Z [J, K] = J J[ @tp c (x) expli J d A x [jSf (x) + J a (x) <p a (x) 



+ - I d x I d y[(p a (x)K ab (x, y) (p b (y)] 
c c 



The connected n-point Green's functions are obtained by taking the n-th func- 
tional derivative of 

W[J,K] = -ilog(Z [J,K]) 

with respect to J. In particular, for n = 1 we obtain the expectation value of the 
field $ a (x): 

However, the first derivative of W [J, K\ with respect to K gives the full 2-point 
Green's function, including disconnected diagrams: 



SW [J, K] 
5K ab (x, y) 



-(G ab (x,y) +<f) a (x)<f) b {y)) . 



[2.1.2) 



J=K=0 



Using Eqs. (j2.1.1|) and ()2.1.2|) . one can express J and K in terms of and G. 
The 2PI effective action V [0, G] is the double Legendre transform of W [J, K] : 



T [0, G] = W [J, K] - J d 4 x [J a (x) a (x)} 

c 



1 



d 4 x / d 4 y 



G ab (x, y) K ab (x, y) 



c c 
+ a (x) K ab (x, y) 4> b (y) 



2. 1 2PI Effective Action 
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From this, we see that 

w [0, gj 

5(f) a (x) 

and 



J a (x) - / d 4 y [K ab (x, y) 6 (y)} 



5T[0,G] 1 ( 

--K ab (x,y) . 



SG ab (x, y) 2 

In the case of vanishing external sources J and K the equations of motion for 
and G read 

5Y [0, G] 



5(f) a (x) 

6T [0, G] 



0, (2.1.3) 
dG afe (x, y) 

The usual 1PI effective action F [(f)] is obtained by Legendre transforming the 
generating functional for connected n-point Green's functions with a vanishing 
bilocal source. Having this in mind, one sees the connection between T [0] and 
r[0,G]: 

r[0] = r[0,G s ^] , (2.1.5) 

where G sta t satisfies the stationarity condition ([2.1.4)) : 



[0, G] 

SG ab (x, y) 



= 

G=G stat 



In order to solve Eqs. (|2. 1.3(1 and (|2. 1.4(1 it is necessary to find an explicit expres- 
sion for T [0, G\. John M. Cornwall, R. Jackiw and E. Tomboulis provided the 
following parametrization in Ref. [3*2] : 

T[0,G] = SM + ^trcflog^jj+^trc^G] (2.1.6) 
+r 2 [0, G] + const . 

Here the squared brackets indicate that the traces, the logarithm and the product 
Gq l G have to be evaluated in the functional sense. S [0] is the classical action 
defined on the closed Schwinger-Keldysh time contour and 

■ r -i ( , 5 2 sm 



5<fi a (x) 5(j) b (y) 
-D x - m 2 - ^4>d (x) 4> d (x) ) 5 ab 5 c (x - y) (2.1.7) 

- 7^0a (x) <p b (x) 5 C (x - y) 
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is the inverse classical propagator. The constant can be chosen such that Eq. 
(I2.1.5|) is indeed satisfied. To compute r 2 we define a new field which has van- 
ishing expectation value: 

K {X) = $a (X) - <Pa {X) • 

Then we rewrite the classical action in terms of $^ (x) and 4> a (x): 

S[$] = S [$' + 0] . 

Apart from linear and quadratic terms, the action on the right hand side includes 
terms cubic and quartic in the field (x). T 2 [</>, G] is the sum of all two-particle 
irreducible vacuum graphs with vertices determined by those cubic and quartic 
terms, while internal lines are set equal to the connected propagator G a b (x,y). 

2.2 1/N- Expansion of the 2PI Effective Action 

In a free field theory r 2 vanishes and one can compute the 2PI effective action 
exactly. However, this is not possible for a theory that contains interactions. 
When this is the case, one has to look for an appropriate approximation scheme. 
In this section, we review briefly the 1/iV-expansion of T 2 to next-to-leading order 
(NLO). Here, we restrict ourselves to the presentation of the results only. Their 
derivation is described in more detail in Refs. ^\ \'27\. 

To obtain the 1 /iV-expansion of T 2 [4>,G] we have to sort all the diagrams 
contributing to T 2 [<p, G] with respect to the powers of their 1/N factors, beginning 
with the lowest powers and proceeding to higher powers. We take all diagrams 
into account that contribute to next-to-leading order, i.e. we write r 2 [4>, G] in 
the form 

T 2 [<j>, G] = if 0) [0, G] + rf LO) [cf),G] + ... (2.2.1) 

and neglect terms beyond next-to- leading order in 1/N. Fig. 12.11 shows the only 
diagram that contributes to leading order. Hence: 

if 0) [0, G] = J d 4 x [G aa (x, x) G bb (x, x)} . 

c 

The next-to-leading order contribution is shown in Fig. 12.21 A bulky vertex indi- 




Figure 2.1: Leading order contribution to T 2 . 



2.2 1/N -Expansion of the 2PI Effective Action 
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Figure 2.2: Next-to-leading order contribution to T 2 - The bulky 
vertices are defined in Fig. 12.31 Crosses indicate an expectation 
value of the original quantum field. 




Figure 2.3: Resummation scheme for the next-to- 
leading order contribution to T 2 [(f), G\. 



cates an infinite series of diagrams, where each diagram has one vertex and one 
loop more than its predecessor. The corresponding recursively defined resumma- 
tion scheme is shown in Fig. 12.31 We obtain 

rf LO) [0,G] = ~tr c [log [5]] 

+ 6^ / ^ X J ^ ^ ^' ^ ^ a ^ Gab ^' ^ ^ b ^ ' 
c c 

Here 

A % A f 

I {x, y) = -^ G ^b (x, y) G ab (x, y) - — J d A z [I (x, z) G ab (z, y) G ab (z, y)] 

c 

(2.2.2) 

resums the chain of bubbles caused by one bulky vertex. For the diagrams that 
contain exclusively four-point vertices this resummation gives the trace over the 
functional series of the logarithm of 

B (x, y) = 5 C (x - y) + ^G ab (x, y) G ab (x, y) . (2.2.3) 

In the symmetric regime (0 = 0) all diagrams that contain three-point vertices 
vanish and only the trace remains as the NLO contribution to T 2 [<fr, G}. 

Eventually, we would like to illuminate the distinction in the resummation 
between the 2PI l/iV-expansion and a perturbative expansion. Altogether, there 
are four diagrams that contribute to T 2 [0, G] to order A 2 , but only two of them 2 

2 These are the second and the third diagram shown in Fig. 12.21 with, according to Fig. 12.31 
the bulky vertex replaced by the classical one. 
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contribute to the 2PI 1/iV-expansion to next-to- leading order. Rotating one of 
the two vertices in each of these two diagrams by 90° respectively, one obtains 
the remaining two diagrams, which contribute at NNLO in 1/N. Similar con- 
siderations hold for higher loop diagrams. At each order in the coupling there 
are only two diagrams contributing to the 2PI 1 /iV-expansion to next-to-leading 
order. 



2.3 Evolution Equations 

From here on, for simplicity we restrict our considerations to the symmetric 
regime, where 

(f) a (x) = . 

The starting point for the derivation of the evolution equations for the spectral 
function and the symmetric propagator is the equation of motion (j2.1.4|) . With 
the parametrization (|2.1.fij) it reads 

C^) = C -^)- 2 yii|L_, (2.3.1) 

which is nothing but the exact Schwinger-Dyson equation for the propagator 
where the proper self energy is given by 

E ab (x,y) = 2i ^ 2 f ] . (2.3.2) 
bG ab (x,y) 

When we insert the inverse classical propagator (|2.1.7|) (here for = 0) as well 
as Eq. (|2.3.2j) the Schwinger-Dyson equation (|2.3.1|) takes the form 

Gab V) = - (°* + ™ 2 ) Sabh (x-y)- Y, ab (x, y) . (2.3.3) 

At next-to-leading order in the 1/iV-expansion of the 2PI effective action Eq. 
()2.3.2|) yields for the proper self energy: 

S«6 (x, y) = -^S ab S c {x - y) G cc (x, x) - j^B' 1 (x, y) G ab (x, y) . (2.3.4) 

By multiplying Eq. (J2.2.3|) from the left with B~ x and using the identity 

J d 4 z [B- 1 (x, z) B (z, y)] =5 c (x- y) (2.3.5) 



we find for B 



-i. 



B- 1 (x, y)=6 c (x-y)-^ J d 4 z [B- 1 (x, z) G ab {z, y) G ab (z, y)] . (2.3.6) 
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a a a b 



Figure 2.4: Diagrammatic expansion of the self energy in the sym- 
metric regime. The bulky vertex represents the resummation that 
has been introduced in Fig. 12.31 Vertices with external lines are 
linked to space time variables which must not be integrated over. 



From this we see, that we can write B l (x, y) in the form 

B- 1 (x, y) =8 c (x-y)- il (x, y) , (2.3.7) 

where I(x,y) is given by Eq. ()2.2.2|) . Further simplification arises if we take into 
account that for systems with 0(iV)-symmetric initial conditions the propagator 
satisfies 

G ab (x,y) = G(x,y)5 ab . (2.3.8) 

Eq. ()2.3.4)1 shows that this property is inherited by the proper self energy, and 
additionally it suggests, together with Eq. ()2.3.7j) . to decompose the self energy 
into a local and a nonlocal part. 

S (x, y) = -?£( local ) ( x ) 5 C (x-y) + £( nonlocal ) ( x , y) 

The diagrammatic expansion of the proper self energy is illustrated in Fig. 12.41 
The two tadpoles in the self energy result from the LO and NLO double bubbles 
in T2 and form the local part of the self energy. Correspondingly, the setting 
sun forms the nonlocal part. While the nonlocal part takes direct scattering 
into account, one immediately sees that the local part causes a mass shift only, 
wherefore we define the effective mass by 

M 2 (x) = m 2 + S (local) (x) 

= rn 2 + \^-^G(x } x) . (2.3.9) 
bJM 

When we replace m 2 in the Schwinger-Dyson equation ()2.3.3)1 by M 2 (x), the 
proper self energy reduces to its nonlocal part which results in 

s( noniocai) (a . )y) = --^G (x,y) I (x, y) , (2.3.10) 
where due to Eqs. ()2.2.2j) and ()2.3.8j) I(x, y) satisfies 

fd 4 z[I(x,z)G(z,y)G(z,y)} . (2.3.11) 



I (x, y) = ^G (x, y) G (x, y) - i^ 
o o 
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Next, we use the analogue to Eq. ()2.3.5|) for the propagator G(x, y) to bring the 
Schwinger-Dyson equation into a form that is appropriate for an initial value 
problem: 

-iSe (x-y)= (□, + M 2 (x)) G{x,y) + iJ d 4 z p nonlocal ) ( x , z ) G {z, y)] . 

c 

(2.3.12) 

Now, we decompose G and £( nonlocal ) according to Eq. (jl.1.3)) . with functions 
that are defined on the closed time path C. 

G (x, y) = 6 C (a; , y°) G> (x, y) + 9 C (y°, x°) G< (x, y) (2.3.13) 



^(nonlocal) ^ y) = ^ y 0) S> y) + Q( , ^ ^ S< ^ y) (2 _3_ 14) 

Using that in the sense of distributions 

□xG (a, y) = C y°) (x, y) + 6 C x°) U X G < (x, y) - iS c (x - y) , 



we find that because of the decompositions ()2.3.13|) and (|2.3.14|) . and the analy- 
ticity properties of G > (x, y) and G< (x, y) with respect to their time arguments, 
Eq. ()2.3.12|) splits in two equations, one for G> (x, y) and one for G< (x, y) 3 : 



r 

(n x + M 2 (x))G > (x,y) = 2jd 3 z J dz°[lm(E > (x,z))G > (z,y)\ (2.3.15) 



and 



dz° [£> (x,z)Im(G > (z >3 /))] 



r v 

(□ ;E + M 2 (a;))G' < (x,y) = 2 J d 3 z J dz° (x, z)lm(G < (z,y))} (2.3.16) 



a; 

- J dz° [Im(S< (x, z))G< (z, 



It is worth noting that in the literature these equations are known as the Kada- 
noff-Baym equations, and that they are used for the description of an enormous 
variety of different phenomena ranging from the very small as in our case over 



3 Here we also use the hermiticity of G> (x,y) and G < (x,y). The same relations hold for 
the proper self energy: E> (x, y) = £> [y, x) = £< (x, y) 
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conductors and semi-conductor devices to the very large as in cosmology. Due 
to the hermiticity property of G>(a;, y) and G < (x,y) one can reduce these two 
complex- valued evolution equations to two real- valued evolution equations. These 
two real-valued evolution equations are the evolution equations for the spectral 
function and the symmetric propagator which we mentioned earlier. Actually, 
the spectral function g(x, y) is a purely imaginary function. Therefore we define 

p(x,y) = ig(x,y) , 

which is obviously a real quantity. The evolution equation for the symmetric 
propagator is given by the sum of Eqs. (j2.3.15|) and (j2.3.16|) . the one for the 
spectral function is obtained from their difference. 



(□ X . + M 2 (x))F(x,y) 



- / d 3 z 



dz°[E p (x,z)F(z,y)] (El) 



dz° [S F (x,z)p(z,y)] 



(n x + M 2 (x)) p (x, y) = J d 3 z j dz° [E p (ar, z) p (z, y)] 

Because of Eq. (12.3. 10(1 and the decomposition ()2.3.14|) 4 

A 



and 



(ar, S/) 
S< (x,y) = 



3iV 
A 



G> (x,y)I> (x,y) 



3N 



G< (x,y)I< (x,y) . 



For the resummation functions I>(x, y) and J<(x,y) we obtain 
A 



I> (x, y) 



6 



G> (x,y) G> (x,y) 



'4 Uh 



(E2) 



(2.3.17) 
(2.3.18) 



dz° [(/> (x, z) - I< {x, z)) G> (z, y) G> (z, ; 



- J dz° [/> (x, z) (G> {z, y) G > (z, y) - G < (z, y) G < (z, y))} 



4 We also decompose the resummation function I (x,y) according to 
I (x, y) = 6 C (x°,y°) /> {x, y) + 9 C (y°, x°) Z< (x, y) 
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and 



A 



I<(x,y) = -G < (x,y)G < (x,y) 
o 



iX 
~6~ 



d 3 z 



dz° [(/> (x, z) - /< (x, z)) G< (z, y) G< (z, y)] 



o 



dz [/< (x, z) (G > (z, y) G > (z, y) - G< (z, y) G < (z, y))} 



The various self energies that appear in the evolution equations (|K1|) and (|E2 
are then given by 

E F (x,y) = - (£> (x, y) + £< (x, y)) 



and 



~ (f (x, y) I F (x, y) - ^p (x, y) I p (x, y) 



E p (x,y) = i (£> (a:, y) - £< (a;, y)) 

= —gft ( F v) h y) + p ( x , y) ^ y) j , 



(E3) 



(E4) 



where the symmetric and the spectral resummation functions If(x, y) and I p (x, y) 
are defined by 

1 



If (x, y) 



(I> (x, y) + I< (x, y)) 
1 



6 



- [F (x,y)- -p (x,y) 



A 
6 









J d 3 z 




I p (x, z) ^ 




L 





and 



dz° [I F (x, z) F (z, y) p {z, y)\ 



I p (x,y) = i(I > (x,y)-I < (x,y)) 
= -^F(x 1 y)p(x } y) 



^Jd'zJ dz°{Ip(x,z)F(z,y)p(z,y)] . 



(E5) 



(E6) 
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The evolution equations (jEl|) to (jE6j) form a closed set of equations which deter- 
mine the evolution of the Schwinger-Keldysh propagator completely. It is worth 
noting that these equations are explicit in time. This is obvious for the self en- 
ergies. For those equations containing integrals it means that, when evaluating 
such an integral for given x° and y°, one only needs to know all the functions in 
the integrand for times between the initial time and x° or y°, respectively. This 
allows for an efficient numerical solution of these equations for any given initial 
conditions. A detailed discussion of this as well as numerically obtained solutions 
for miscellaneous initial conditions can be found in Ref. 

We mentioned in the introduction that the 1/iV-expansion of the 2PI effective 
action is capable of controlling nonperturbatively large fluctuations. Now, we 
are able to prove this statement PQ I2|- in the case of nonperturbatively large 
fluctuations 



One sees immediately, that this property is inherited by the resummation func- 
tions If (x, y) and I p (x, y). As a consequence every diagram contributing to the 
1/iV-expansion of the 2PI effective action to next-to- leading order in the symmet- 
ric regime contributes with the same order in the coupling. Therefore any loop 
expansion breaks down in such a situation and one has to use a nonperturbative 
approach, as provided by the 1/A-expansion of the 2PI effective action. 

In the next chapter, using a combination of a first order gradient expansion 
and a Wigner transformation, we are going to derive kinetic equations from these 
evolution equations. In particular, by additionally employing a quasi-particle 
(on-shell) approximation, we are going to show how one can obtain a generalized 
Boltzmann equation. In this sense one can consider the full evolution equations 
as quantum Boltzmann equations resumming the gradient expansion up to infinite 
order and including off-shell effects [33]. 




and 
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Chapter 3 

Local Approach to Equilibrium 



In this chapter we consider quantum fields in the intermediate-time drifting and 
the late-time thermalization regimes. As pointed out in the introduction we will 
exploit the slight dependence of the propagator on the center coordinates and 
the effective memory loss. By combining a first order gradient expansion and a 
Wigner transformation we derive kinetic equations that are local in time from the 
full nonlocal evolution equations. These kinetic equations are important tools for 
the description of the evolution of nonequilibrium quantum fields for large times 
where the memory integrals complicate computations. 

Additionally employing a quasi-particle approximation, we obtain a general- 
ized Boltzmann equation. In the derivation of the classical Boltzmann equation 
one assumes that the system is so dilute such that there are practically only colli- 
sions which involve two particles and that one can neglect collisions between three 
or more particles. In contrast to the classical Boltzmann equation our generalized 
Boltzmann equation includes an effective mass and an effective coupling which 
arise from the resummation of the infinite series of diagrams that contribute to 
the 1/iV-expansion of the 2PI effective action to next-to-leading order. Due to 
the effective coupling we expect our generalized Boltzmann equation to be valid 
even for systems which are so dense that one cannot neglect collisions between 
three or more particles. 



3.1 Quantum Kinetic Equations 

We begin the derivation of the kinetic equations with the evolution equations for 
the symmetric propagator and the spectral function (jEl|) and (|E2ji . namely 




(3.1.1) 



/ 



d 4 z 9 (z°) (p R (x, z) F (z, y) + S F (x, z) G A (z, y) ) 
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and 



[o x + M 2 (x))g(x,y) 



(3.1.2) 



d 4 z 



9 (z°) (x, z) G A (z, y) + E R (x, z) g (z, y) j 



Now, we interchange x and y in Eq. ()3.1.1|) and use the symmetry properties of 
the symmetric propagator and the self energies with respect to the order of their 
arguments. Subtracting the obtained equation from Eq. ([3.1.10 . we get: 

(n x - □„ + M 2 (x) - M 2 (y))F(x, y) (3.1.3) 

= J d 4 z 6 (z°) (f (x, z) S a (z, y) + G R (x, z) S f (z, y) 

- Y< R (x, z) F (z, y) - S F (x, z) G A (z, y) 



The same procedure yields for Eq. (j3.1.2|) : 

(n x -n y + M 2 (x) - M 2 (y) ) q (x, y) 



(3.1.4) 



(z°) [G R (x, z) £ e (z, y) + q (x, z) Y, A (z, y) 



- H e (x, z) G A (z, y) - S H (x, z) g (z, y) \ 



For every two points v and w in space time, their center and relative coordinates 
are denoted by 

v + w 



and 



v — w . 



However, the space time coordinates x and y on the left hand sides of Eqs. ([3.1. 3j) 
and (j3.1.4[) are particularly important: 



X = X xy = — — and 



s = s 



■ni 



x - y . 



(3.1.5) 



Every function / of v and w can be written as a function of the center coordinate 
X vw and the relative coordinate s vw : 



f(v,w) = f[ x vw + x vw - 



We now rewrite Eqs. ([3.1.30 and (j3.1.4[) using center and relative coordinates. 
An application of the chain rule gives 



u x - u v = d x „d Xu - dy^dy^ = 2d Sa d_ 
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In the drifting regime, the exponential suppression of correlations with the initial 
time allows us to send the initial time to — oo, i.e. from here on we drop the 
9 function in the above memory integrals. Hence, the expressions on the right 
hand sides of Eqs. (J3.1.3|) and (J3.1.4j) include integrals of the form 



f(x,z)g(z,y) 



d 4 z 



f(X x 



zy: s zy) 



(3.1.6) 



In equilibrium the propagator and the proper self energy are invariant under 
space-time translations. This means that they depend on the relative coordinates 
only, and that there is no dependence on the center coordinates. This does not 
hold for a system out of thermal equilibrium, where the propagator and the 
self energy depend on the center coordinates. However, if the propagator is a 
sufficiently smooth function of the center coordinates — which it is by definition 
in the intermediate-time drifting regime (see Fig. ^) — a Taylor expansion to first 
order in the center coordinates should be a good approximation. At this point 
we would like to emphasize that the system still may be far away from thermal 
equilibrium and that we do not justify the Taylor expansion by assuming that the 
system has approached the equilibrium sufficiently closely. The Taylor expansion 
of the mass terms to first order in dx? around A gives 



M 



x+ s - 

2 



M 



x- s - 

2 



s» ■ (d x ,M 2 {X)) . 



In the integral (J3.1.6|) the expansion of / and g into a Taylor series to first order 
in dx^ around X xz = X and X zy = X respectively, yields 



d A z 



f(x,z)g(z,y) 



d A z 



f (X, s xz ) g (X, s zy ) - -ff (A, s xz ) (d x ,g (X, s zy )) 



+ -y [dxrf (A, s xz )j g (X, s zy ) 

After we have removed the 9 functions from the memory integrals and performed 
the first order gradient expansion in Eqs. (J3.1.3|) and (J3.1.4j) . our next step is to 
Fourier transform these equations with respect to s. On the left hand side we get 



LHS&LQ 

d 4 s 



exp (iks) (2d Sii d x » + s M (d x ,M 2 (X))) F (A, s) 



-i (2Wd x » + {d x ,M 2 (A)) d h ) F (A, k) 



and similar for Eq. (|3. 1.4)1 . F(X,k) denotes the so-called Wigner transform of 
F (x, y), which is nothing but the Fourier transform of F (A, s) with respect to 
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s. It is worth noting that due to the hermiticity of G> (x,y) and G< (x,y), the 
Wigner transforms of the symmetric propagator and the spectral function are 
real functions. 

On the right hand side of Eqs. (|3.1.3|) and (|3.1.4|) . we find two different types 
of integrals, which we transform individually: 



d A z 



f(X, s xz )g(X,s 



zy) 



d 4 s / d 4 z 



exp (iks) f (X, s xz ) g (X, s zy ) 



f(X,k)~g(X,k) 



and 



d A z 



-^f (X, s xz ) (d x ,~g (X, s zy )) + -^(d x >f (X, s xz ) ) g (X, s zy} 



d s Id z 



exp (iks) [ -^f(d x ,~g) + 



- ((dkj (X, k)j (d x ,~g (X, k)) - [d x ,f (X, k)j (d k J (X, k) 



f(x,k)-~g(x,k) 



PB 



In the last step, we have introduced the Poisson brackets to shorten the notation. 
After all these transformations Eqs. (j3.1.3|) and (|3.1.4j) become 

-i (2k»d x » + (d x ,M 2 (X)) d kfl ) F (X, k) 



(3.1.7) 



-FIE 



-*R 



+ E F G 



<R 



G t 



+ - { F; E R + E A 



— — i Ep\ Gp> 
pb 2 



G, 



PB 



and 



(2k»d x „ + (d x ,M 2 (X)) d K ) q (X, k) 



(3.1.8) 



E n [G 



G, 



Q [E 



-^lE e ;G 



' R 



G, 



J R 



PB + 2 



R 



PB 



where all the functions on the right hand side depend on X and k. From Fourier's 
theory, we know that for two functions / (X, s) and g (X, s) the Fourier transform 
of their product with respect to s is given by the convolution of their individual 
Fourier transforms 

d A p 



d*s 



exp (iks) f (X, s) g (X, s) 



(2vr) 4 

f*§) (x,k) 



f(X,k-p)g(X,p) 
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and that the Fourier transform of a 6 function in one dimension is given by 



dx [exp (ikx) 9 (x)} = lim 



o k + ie 

Hence, we can write the Wigner transforms of the retarded and the advanced 
propagator in the following form: 



(3.1.9) 
(3.1.10) 



r dp 


p(X,k,p°)l 


J 27 


k° — p° + ie 


r dp 


lg(X,k,P°)] 


J 2^ 


k° — p° — ie 



As the 5 function can be approximated by 



5e (k°) 



we find that 



7T {kl + #) ' 

G R (X, k)-G A (X, k)=ig(X, k) 



(3.1.11) 



The same relation holds for the proper self energy. On the other hand, we have 
already seen that g (X, k) is a real function. Therefore, Eqs. ()3.1.9|) and (J3.1.1()j) 
yield 

G* R (X,k)=G A (X,k) . 



From this we see immediately that 

G R (X, k) + G A (X, k) = 2 • Re (g r (X, k) 

and 

g (X, k) = 2 ■ Im (Gr (X, k) 



(3.1.12) 
(3.1.13) 



Again, the same relations hold for the self energy, and Eqs. (|3. 1.7(1 and ()3.1.8() 
now have been simplified to take the following form: 



(2Pc>xm + {d x .M 2 (X)) d fc J F (X, k) 

= F (X, k) t Q (X, k) - g (X, k) t F (X, k) 



(Kl) 



+ \ S F (X, k) , Re ( G R (X, fe) ) + | Re (^E R (X, fc)J , F (X, fc) j ^ 



and 



(2A;^ X m + (9 X mM 2 (X)) 9 fc J g (X, k) (K2) 
X (X, fc) , Re (Gr (X, |Re {t R (X, fe)) , g (X, fe)' 
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These are the kinetic equations for the symmetric propagator and the spectral 
function which we mentioned earlier. In the next section we will show that, 
under certain assumptions, one can derive a generalized Boltzmann equation 
from Eq. (jKlj) . As we will see, the left hand side of Eq. (jKl|) will reduce to 
the usual form of the left hand side of the classical Boltzmann equation as for 
example presented in the third section of Ref. [22], and on the right hand side 
FT, e — qYif will evoke a collision integral. 

The convolution of the Wigner transformed spectral function with the Wigner 
transformed 9 function as seen in Eq. (J3.1.9j) leads to a principal value integral 
which in general cannot be evaluated. Therefore, we cannot use Eq. (J3.1.9J) to 
calculate the Wigner transform of the retarded propagator, which appears in the 
kinetic equations (|K1J) and (jK2|) . We circumvent this difficulty by deriving an 
extra equation for the Wigner transform of the retarded propagator. For the same 
reason we will derive a similar equation for the Wigner transformed retarded self 
energy as well. We obtain the equation for the Wigner transform of the retarded 
propagator in the following way: By calculating the effect of O x on Gr (x, y) and 
Ga (x, y) we find 

i6 (x° - y°) n x g (x, y) = n x G R (x, y) - S (x - y) , 

-id (y° - x°) D x g (x, y) = a x G A (x, y) - 5 (y - x) . 

From this we see that the multiplication of Eq. (j3.1.2j) once with i9 (x° — y°) and 
once with —id (y° — x°) yields equations for Gr (x, y) and Ga (x, y) respectively: 

5(x-y)= (D x + M 2 (x)) G R (x, y) + J d^z [T. R (x, z) Gr (z, y)\ (3.1.14) 

and 

6{y-x) = {a x + M 2 (x))G A (x,y) + J d 4 z[Z A (x, z) G A (z,y)} . (3.1.15) 

We add Eq. (J3.1.15j) . with x and y interchanged, to Eq. (J3.1.14|) . which gives 

28{x-y) = (a x + O y + M 2 (x) + M 2 (y))GR(x,y) (3.1.16) 
+ J d 4 z[T l R(x,z)G R (z,y) + G R {x, z)Y, R (z,y)} , 

and proceed along the lines that lead from Eqs. (|3.1.3|) and (|3.1.4p to the kinetic 
equations ()K1|) and (jK2|) . When, we rewrite Eq. ()3.1.16|) using center and relative 
coordinates and perform a first order gradient expansion, we find that the effective 
masses give 

M 2 (x) + M 2 (y) = 2M 2 (X) . 
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Furthermore, neglecting terms beyond first order in dx^, we obtain 



2d s nd s 



The integral on the right hand side behaves exactly in the same way as the ones 
that we met during the derivation of the kinetic equations (jKlj) and (jK2|) . Finally 
we perform a Fourier transformation with respect to s. The result is an algebraic 
equation for Gr (X, k): 

-k 2 + M 2 (X) + t R (X, fc)) Gr (X, k) . (K3) 

Due to Eqs. (jK3|) and ()3.1.11jl . we can express the spectral function in terms of 
the proper self energy: 



g(X,k) 



(X,k) 



-k 2 + M 2 (X) + Re (E R (X, A;)))' + (\t e (X, k) * 



(3.1.17) 



We still need to work out the expressions of the self energies that appear in the 
kinetic equations. In addition to the expressions (jE3|) and (jE4|) for the symmetric 
and the spectral self energy we find for the retarded self energy: 



Zr (x,y) 



id {x°-y°) E e (x,y) 
A 



3X 



(F (x, y) I R (x, y) + Gr (x, y) I F (x, y) ) . (3.1.18) 



The Wigner transformation of Eqs. (|E3|) . (|E4|) and ()3.1.18|) yields: 



S F (X, k)~- 

t Q {X,k) 
t R (X, k)-- 



X 



3N 
A 

~3X 

^_ 

"3X 



F*I F )(X,k) + -[Q*I e )(X,k) 



F*I e ) (X,k)+(g*I F ) (X,k) 



F*lR)(X,k)+{G R *I F )(X,k) 



(K4) 

(K5) 
(K6) 



In the derivation of the kinetic equations we exploited the effective memory loss 
in order to drop the 9 functions that appear in the full evolution equations (|3.1.1j) 
and ()3.1.2j) . By the same reasoning we can drop the respective 6 functions in the 
equations for I F (x,y) and I g (x,y) which then read: 

If{x,d) = ^ (f 2 (x,y) + ^e 2 (x,y) 



al A z 



/„(./■.:) ( F 2 (z,y) + ^g 2 (z,y) 



+ 2I F (x,z)F(z,y) G A (z,y) 



(3.1.19) 
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and 



i e ( x >y) = 3 F ( x > v) q ( x > v) 



I Q (x,z)F(z,y) G A (z,y) 



+ I R (x,z)F(z,y)g(z,y) 



(3.1.20) 



Multiplying Eq. ()3.1.20|) with id (x° — y°) gives an equation for the retarded re- 
summation function I R (x, y): 

I R (x,y) = ^G R (x,y)F(x,y) (3.1.21) 



J d A z [I R (x,z)G R (z,y)F{z,y)} 



When we rewrite Eqs. fl3.1.19j) . ()3.1.2()j) and ()3.1.21|1 using center and relative 
coordinates, expand the functions inside the integrals to first order in dx» around 
X and Fourier transform the equations with respect to s, we obtain the following 
equations for the Wigner transformed resummation functions: 



i F (x,k) = ^(7f*f) + 



(K7) 



I R - ( [F*F] +-(g*g) ) +2I F [F*G A 



+ l{f R -\F*F) + 1 -(^~ g) ^ B + i { I F ^F*G A/jpB 



Ir (f*q)+ I e (P * G A ) (K8) 

I R (x,k) = ±(G R *F s ) ( 1 -Ir) + j{(gr*f);Ir} pb ■ ( K9 ) 

Of course, the functions and convolutions on the right hand sides depend on the 
center coordinate X and the momentum k, without exception. 

We stress that the kinetic equations ()K1|) to (jK9|) form — exactly as the evolu- 
tion equations ()E1|) to (jE6|) in Sect. 12.31 — a closed set of equations. Their range 



= 3(^*^-3 
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of applicability is restricted only by the gradient expansion that we employed 
in their derivation. Hence, in regimes where this gradient expansion is justi- 
fied, the kinetic equations (jKlj) to (jK9|) describe the evolution of nonequilibrium 
quantum fields without further approximations. In this context we would like to 
emphasize that our kinetic equations indeed inherit the capability of describing 
nonperturbatively large fluctuations from the nonlocal evolution equations, as it 
was described in Sect. 12.31 j2J. Furthermore, it is important to note that we de- 
rived our kinetic equations starting from the effective action. As a consequence 
it can be shown that our kinetic equations possess the generic feature of exact 
conservation laws at the level of expectation values |3"E] . 

Before we close this section, we convince ourselves that our kinetic equations 
contain a consistent description of the equilibrium. We mentioned already in 
Sect. I1.2l that in equilibrium the propagator, and thus the proper self energy too, 
are space-time translation invariant. Hence, they do not depend on the center 
coordinate X, and so the left hand sides of Eqs. (jKl|) and (jK2|) as well as the 
Poisson brackets on their right hand sides vanish. From this we see already that 
at least Eq. (jK2|) is identically fulfilled in equilibrium. To see that Eq. (|K1|) also 
describes the equilibrium in a consistent way, it remains to show that the collision 
term also vanishes in equilibrium. To see this most easily, we note that we can 
write the collision term also as 

F (k) t Q (k) - q (k) t F (k) = G< (k) £> (k) - G> (Jfe) £< (k) . (3.1.22) 

For the same reason as above, the Poisson brackets in Eqs. (jK7|) and ()K8|) vanish, 
and these equations can be brought into the following form: 



= £ (G> * G>) (*) 1 7 ! R (fc \ ■ (3. 1 -23) 

V / i I A 

1 " r 3 



(G*R * F) (k) 



I< (*) = \ (G< * G<) (k) • *,* R{k \ • (3-1.24) 

6V 1 l + l(G* R *F)(k) 

With the aid of these equations we can show that in equilibrium the proper self 
energy satisfies 

£< (k) = exp (-pk°) S> (k) , 

and together with the Fourier transformed Kubo-Martin-Schwinger condition 
(|1.2.9|) for the propagator we see then immediately, that in equilibrium the colli- 
sion term in the kinetic equation for the symmetric propagator indeed vanishes 
and that the kinetic equations really describe the equilibrium consistently. 



3.2 Generalized Boltzmann Equation 

The first step to learn something about the kinetic equations that we derived 
in the last section is to consider the special case of a system containing long- 
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living single-particle excitations, where the width of the spectral function is much 
smaller than the effective particle mass. This allows us to employ a so-called 
quasi-particle approximation, and we can derive a generalized Boltzmann equa- 
tion from the kinetic equations (jKlj) to (jK9|) . Our starting point is the so-called 
Kadanoff-Baym ansatz, which is suggested by Eqs. ()1.2.10|) and ([1.2.110 . namely: 

G< (X, k) = g (X, k) n (X, k) , 
G> (X, k) = g (X, k)(l+n (X, k) 



It follows then that 



F(X.k) n(X.k) [n(X,k) + ^] . (3.2/1) 



To employ a quasi-particle approximation means to assume that the particle 
number density does not depend explicitly on the energy, i.e. 

d k on (X, k) = , (3.2.2) 

and that the spectral function g (X, k) has the same form as for a free theory: 

g (X, k) = 2vr • sign (A; ) • 5 (k° 2 - E 2 (X, k)) 

where 

E (X, k) = Vk 2 + M 2 (X) . 

We see that we obtain this form of the spectral function from Eq. (J3.1.17[) when 
we set T (X, k) = — ^Ti e (X, k) and consider the limit Y (X, k) — > 0. However, we 
have to emphasize the fact that we cannot send S e (X, k) — > everywhere, as 
for example in the collision term, where the spectral self energy is necessary for 
thermalization. One can verify directly that in the sense of distributions 

(2k»d x » + (d x ,M 2 (X)) d kl ) q (X, k) = . (3.2.3) 

An immediate consequence from Eq. (j3.2.3|) is the fact that we can employ the 
quasi-particle approximation consistently only, if the Poisson brackets on the right 
hand side of Eq. (jK2|) vanish. The terms that contain derivatives with respect 
to k° vanish according to Eq. ([3.2.2)1 . and later we will consider spatially homo- 
geneous systems, where also the terms with spatial derivatives vanish. However, 
for the moment we just simply assume that all Poisson brackets vanish, although 
we consider a spatially inhomogeneous system. Inserting Eqs. (13.2.1)) and ()3.2.3j) . 
the left hand side of Eq. ()K1)) becomes 

LHSm = E k) (S (k° - E (X, k)) - 6 (k° + E (X, k))) (3.2.4) 

x (2k?dx» + (d x ,M 2 (X)) <9 fe J n (X, k) . 
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As we neglect the Poisson brackets, due to Eq. (j3.1.22|) and the Kadanoff-Baym 
ansatz the collision term on the right hand side reads: 

RHSm = E { l k) (S (k° - E (X, k)) - 6 (k° + E (X, k)) ) (3.2.5) 
x (£> (X, k) n (X, k) - £< (X, k) (1 + n (X, fc))) . 

On the mass shell we define 

n(X,k) = n (X, k, E (X, k)) , 

which can be interpreted as a phase space distribution function for quasi-particles 
with momentum k and energy E (X, k). Equating the positive energy components 
of Eqs. ()3.2.4|) and (j3.2.5j) and integrating over k° gives 1 : 

(*" + £<tk) s * + wkv {dx ' Ml {x)) dt ) n (x ' k) 

= 2E L k) (S> (X, k) n (X, k) - S< (X, k) (1 + n (X, k)) ) . (3.2.6) 
When we use the notation 



E(X,k) 
and the fact that 

1 



2E (X, k) 



(d X3 M 2 (X)) d k . = -V 1 (d X j E (X, k)) d, 



we can simplify the left hand side. On the right hand side we use the Wigner 
transforms of Eqs. (|2.3.17j) and (|2.3.18j) together with the identities 



n(X,-k) = -l-n(X, k) 
J d 3 p [f (p)] = J d 3 p [f (-p) 



1 This is the step where it is necessary that the Poisson brackets in the kinetic equation for 
the symmetric propagator vanish. If this were not the case, we could not isolate the positive 
energy components on the right hand side. 
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such that we find a preliminary version of the generalized Boltzmann equation: 

(d x o + v ■ V x - (V x £ (X, k)) ■ V k )n (X, k) (3.2.7) 
A f d 3 p 



i2X£(X,k)y (2 



^7T 



£(X,p) 



x |/> (X, k - p, E {X, k) - E (X, p)) (1 + n (X, p)) n (X, k) 

+ /> (X, k - p, £ (X, k) + E (X, p)) n (X, p) n (X, k) 

- /< (X, k - p, £ (X, k) - £ (X, p)) n (X, p) (1 + n (X, k)) 

- /< (X, k - p, E (X, k) + E (X, p)) (1 + n (X, p)) (1 + n (X, k)) 

This equation resembles very much the classical Boltzmann equation as it is 
presented for example in Ref. |22j . However, there are some subtle differences 
that will be explained later. The comparison between the left hand side of this 
equation and Eq. (3.3) in Ref. |22 a shows that E (X, k) acts as a potential for the 
quasi-particles. 

Now, after having pointed out the relation between the left hand sides of 
the kinetic equation for the symmetric propagator and the classical Boltzmann 
equation, we specialize to spatially homogeneous systems, i.e. all quantities lose 
their dependence on the spatial center coordinate X. As pointed out above, an 
immediate consequence of this is the fact that in this case all Poisson brackets 
vanish automatically. To simplify notation, we write t instead of X° from now 
on. The resummation functions in Eq. ()3.2.7|) are related to the ones we used in 
the last section by 



1 s 
2 
and 

1 ~ 
2 

Hence, we obtain from Eqs. (jK7|) and (|K8|) 



/> (t, k) = i F {t, k) + -i e (t, k) 



J< (t, k) = I F {t, k) - -I e (t, k) 



7> (t, k) = ± (G > * G>) (t, k) 1 lR (3.2.8) 

6V 7 1 + f [G* R *F)(t,k) 

and 

J< (t, k) = ± (G < * G K ) (t, k) • ( 3 - 2 - 9 ) 

6V 7 1 + | ( G *r * F ) (*, k) 

We see that the resummation functions are given by the product of their contri- 
butions to the corresponding self energies at two-loop level (setting sun, which 
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means three-loop level for Y2) times some factor which is the same for both re- 
summation functions. This common factor includes the corrections coming from 
the resummation shown in Fig. 12.31 Hence, we can introduce an effective coupling 
associated with the bulky vertex shown in Fig. 12.41 



A e ff (t, k) 



X[l-I R (t,k) 



(3.2.10) 



We will evaluate this effective coupling in the next section. The contributions of 
the resummation functions to the two-loop self energies are given by 



(G>*G>)(f,A0 

1 r d 3 p 

2 J J2nf 



(3.2.11) 



d 3 q 



*(q-(k-p)) 



7T 



E(t,q)E(t,p) 



x <^ 5 (k° - (E (t, p) + E (t, q))) (1 + n (t, q)) (1 + n (t, p)) 



+ 5 (k° 
+ 8(k° 



(-E (t, p) + E (t, q))) (1 + n (t, q)) n (t, p) 
(E (t, p) - E (t, q))) n (t, q) (1 + n (t, p)) 



+ 6 (k° - (-E (t, p) - E (t, q))) n (t, q) n (t, p) 



and 



[G < *G < ){t,k) 

1 r d 3 p 

2 y (2^7 



(3.2.12) 



d 3 g 



*(q-(k-p)) 



7T 



£(*,q)£(t,p) 



X 



5 (fc° - (E (t, p)+E (t, q))) n (t, q) n (t, p) 



+ 5 ( fc o _ ( _£ (t) p) + ^ (t) q)) ) „ (t> q ) (1 + n ( t; p )) 
+ 5 (A; - (E (t, p) — E (t, q))) (1+n (t, q)) n (t, p) 

+ 5 ( fc o _ (_ E ( t) p ) _ £ ( t) q)) ) (1 + n (t; q)) (1 + n ( t) p)) I 



When we insert Eqs. (I3~2~5J) to (13.2. 12jl into Eq. (I3~2~7|) . we obtain on the right 
hand side of Eq. (|3.2.7|) an integral over p and q of a sum of products which 
contain the energy 8 functions appearing in Eqs. (|3.2.11)) and (|3.2.12j) . the effec- 
tive coupling appearing in Eqs. ()3.2.8|) and (|3.2.9|h and the n's and (1 +n)'s of 
Eqs. (j3.2.7jl . ()3.2.11|) and ()3.2.12|) . These products describe the various processes 
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Figure 3.1: Scattering of two particles. 



that can occur, n represents an incoming particle and (1 + n) an outgoing one. 
Depending on the composition of the products, one can classify them to belong 
either to the gain terms or to the loss terms. For example the terms 

n {t, k) n (t, p) (1 + n (t, q)) (1 + n (t, k - p - q)) (3.2.13) 

and 

(1 + n (t, k)) (1 + n (t, p)) n (t, q) n (t, k - p - q) (3.2.14) 

describe the scattering of two particles. The term ()3.2.13|) describes the scattering 
of a particle from the mode k at some particle from the mode p into some other 
mode. Therefore, and because of its emerging with a minus sign, it belongs to 
the loss terms. On the other hand, the term (|3.2.14|) describes the scattering of a 
particle from an arbitrary mode into the mode k. Coming with a plus sign, this 
term belongs to the gain terms. Both processes are illustrated in Fig. 13. if . 

Similarly, the following two terms describe — and Fig. 13.21 illustrates — the 
decay of one particle into three particles and the recombination of three particles 
to one particle, respectively: 

n (t, k) (1 + n (t, p)) (1 + n (t, q)) (1 + n {t, k - p - q)) (3.2.15) 

and 

n(£,k)7i(*,p)7i(t,q)(l + n(t,k-p-q)) . (3.2.16) 

However, momentum conservation prevents the arguments of the respective en- 
ergy 5 functions from being zero. Hence, these decay and recombination processes 
are forbidden and we can cancel these terms in our equation. We prove this state- 
ment for one special case. The energy satisfies the following inequalities: 

E (t, k + p) < ^4M 2 (t) + (k + p) 2 < E (t, k) + E (t, p) . 

The first inequality results from the fact that in the symmetric regime the effec- 
tive mass is always strictly greater than zero and the second one is the triangle 

2 If one finds the arrangement of the momenta inappropriate, because it contradicts the 
natural intuition of momentum conservation, one can use that n and E depend only on the 
magnitude of their momentum arguments and perform a transformation for the momenta, to 
obtain the natural arrangement. 
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Figure 3.2: Decay and recombination processes. 



inequality for the Euclidean norm in M 4 applied to the vectors (M (t) , k) and 
(M (t) ,p). The replacement k — > k — p reveals the following complementary 
inequality: 

E (t, k)-E (t, p) < E (t, k - p) . 
An appropriate choice of the momenta then shows that 

E {t, k) - E (t, p) < E (t, k - p) < E (t, q) + E (t, k - p - q) 

E (t, k) - E (t, p) - E (t, q) - E (t, k - p - q) < 

Hence, the term with the 5 function that has the left hand side of the last in- 
equality as its argument vanishes. With very much the same reasoning one can 
rule out all the remaining decay and recombination terms. 

Furthermore, it is obvious that terms which describe four incoming parti- 
cles leaving to nowhere or four outgoing particles coming from nowhere vanish, 
because the arguments of the respective energy 5 functions have no zeros. 

Eventually, we obtain for the generalized Boltzmann equation 3 : 



d t n (t, k) 



Xn f d 3 p f d 3 q 



1UN 



EkE p -Eq-Ek-p-q 



frr)*J (2tt) 

x (2A efr (t, k - q, E h - E^j + A eff (t, k-p,E k + E t 



c-p-q 



X 



(n k n p (1 + n q ) (1 + n k _ p _ q ) - (1 + n k ) (1 + n p ) r2 q n k _ p _ q j 



In contrast to the classical Boltzmann equation our generalized Boltzmann equa- 
tion includes an effective mass and an effective coupling. The effective mass, 
which in Eq. flBl|) is hidden in the energies, arises from the local tadpole con- 
tributions to the self energy (cf. Fig. 12.4)1 and satisfies a gap equation as will be 
shown in the next section. The effective coupling arises from the resummation of 
the infinite series of diagrams that contribute to the nonlocal part of the self en- 
ergy (cf. Figs. I2.4l and l2.3j) . As we will show in the next section this resummation 

3 We use the following abbreviations: Ey = E (t, k) and = n (t, k) 
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takes higher products of n's and (1 +n)'s into account. Therefore, in contrast to 
the classical Boltzmann equation, we expect our generalized Boltzmann equation 
to be valid even for large particle number densities, where one cannot neglect 
collisions between three or more particles. 

Before we finish this section, we illuminate briefly our expectations about the 
solution of Eq. (JBTJ. As our system contains only bosonic quasi-particles, the 
equilibrium particle number density is the Bose-Einstein distribution function 



n B (k) 



exp (^1) - 1 



(3.2.17) 



In equilibrium there is no time-dependence at all, and 



E(k) 



M% + k^ 



The temperature T of the system is defined only in equilibrium and is determined 
by fitting the numerical data that is obtained at the end of the numerical evolution 
of the Boltzmann equation to the Bose-Einstein distribution. In equilibrium 
the effective mass is a constant and is denoted by M eq . Furthermore, after the 
system has equilibrated, the particle number density must not change anymore. 
Consequently, the collision integral on the right hand side of the Boltzmann 
equation must vanish. We show that this is indeed the case by simply inserting 
the Bose-Einstein distribution function 1)3.2.17)1 into the collision integral on the 
right hand side of the generalized Boltzmann equation (jBlj) 4 : 



d 3 p f d 3 q 



(2tt) 



5 ( E k + E p — Eq — -Ek-p-q 
-^k-^p-E'q-E'k-p-q 



X 



2A cfr (t, k - q, E k - £ q ) + A cff (t, k-p,E k + E p 

x | exp (f3Eq) exp (/3£ k - P - q ) - exp (/3£ k ) exp ((3E P ) 



x 



d 3 p f d 3 q 



(27T) d J (27T) 



5(£ k + £p-£ q -£ 



k-p-q 



EkEpEqEk-p- 



2X cS (t, k - q, E k - £ q ) + A cff (t, k-p,E k + E p 

(e/Jik _ i) ( e ^P - l) ( e ^ q _ [) ( e /3£; k -p- q _ fj 



4 As in equilibrium there is no time dependence at all, here we use the abbreviation E k = 
E (k). is as usual the inverse temperature. 
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x 



( exp ((3 (E k + E P - £ k _ p _ q )) exp (/?£ k - p - q ) - exp ((3 (E k + E p )) 







As was expected, our generalized Boltzmann equation contains a consistent de- 
scription of the equilibrium. In the next section we will derive equations that 
determine the effective mass and the effective coupling for any given particle 
number density at some time t. 



3.3 Effective Mass and Effective Coupling 

Effective Mass 

According to Eq. (|2.3.9|) the effective mass is given by 

M 2 (X) = m 2 + A^^G(X,X) . 

We use that 



G (X, X) — F (X, X) — F (X, 0) 

and express the symmetric propagator in terms of its Fourier transform. The 
quasi-particle approximation then leads to a gap equation for the effective mass, 
which has to be solved self-consistently. 



M 2 (X) 



m 2 + A 



N + 2 f d 4 p 



QN 



2 . ,N + 2 



m z + X- 



12N 



(2vr) 4 
d 3 p 



F{X,p) 

' 2n (X, p) + I 
VM 2 (X) + p 2 



(B2) 



Of course, for a spatially homogeneous system the effective mass loses its depen- 
dence on X, and M 2 (X) becomes M 2 (t). 

Effective Coupling 

Taking into account that Eq. (|K9J) yields 

A (Gr * F) (t, k) 



Ir (t, k) 



we obtain for the effective coupling (|3.2.f (Jj) 

A e ff (t, k) = 



i + Ug r *f) (t,k) 



(B3) 



A 



Gr * F) (t, k) 



(B4) 
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The convolution in the denominator yields 



G R *F)(X, k) 



d 3 p 



2 J (2 



7T 



n(t,p) 



E(t,p) 



G R (t,k-p,k°-E(t,p)) 



G R (t,k-p,k° + E(t,p)) 



where the retarded propagator is given by Eq. ()K3jl as 

1 



G R (t, k) 



-k 2 + M 2 (t) + Z R (t, k) 
Due to Eqs. (|K6|) and (|K7J) the expression for the retarded self energy reads 



(B5) 



t R (t,k) 



A 



d 3 p 



QN J (2 



in) 



n (t, p) + \ 
E(t,p) 



J fl (*,k-p, fc°-.E(*,p)) 



+ 4(t,k-p,A; + £(t,p)) 



(B6) 



A 



d 3 p f d 3 q 



72N J [2-kY J (2 



l-K] 



E(t,p)E(t,q) 



G R (t, k - p - q, k° — E p — £ q ) A eff (t, p + q, £ p + £ q 



n(t,p) + ^J (n(t,q) + ~^ + l - 
+ 2G# (t, k - p - q, k° - E p + £ q ) A cfr (t, p + q, E p - E q 

x f- (n {t, p) + ~J (n (t, q) + + ~ 
+ G fl (t, k - p - q, k° + E p + £ q ) A cfr (t, p + q, -£ p - £ q ) 

n(*,p) + ^j (™(*,q) + ~) +^ 



We find that Eqs. (jBlj) to (|B6j) form a closed set of equations which determine 
the evolution of the particle number density completely. In the next section we 
are going to present methods which allow for the tracing of this evolution for a 
given initial particle number density. 



3.4 Numerical Implementation 

In order to solve the kinetic equations for the symmetric propagator and the 
spectral function or the generalized Boltzmann equation numerically, we propose 
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to use a standard lattice discretization [37]. We consider our whole system of 
physical particles to live on a lattice enclosed in a spatial cube with periodic 
boundary conditions. When doing so, the cube and the periodic boundary con- 
ditions cause the momenta to be discrete. Conversely, because of the lattice the 
momenta are periodic, which means that the lattice takes care of the regulariza- 
tion. Physically reasonable results are obtained by approaching the continuum 
as well as the infinite volume limit. 

The edges of the cube have length L, such that its volume is V — L 3 . The 
lattice spacing is a s , and N s is the number of lattice bins along each edge, such 
that 

L = N s a s . 

For a spatially homogeneous system the lattice discretization indicated above 
yields for the momenta: 

9 ^— \ 4 . r, / a s k n 
k -> > snT 

where as usual, the momenta k n are discretized according to 

k n = ^n, 

with 

nG {0,...,iV s -l} 3 . 

As a consequence, we also have to reverse the thermodynamic limit and perform 
the following replacement: 

d 3 p 



1 E 

v ^ 



(2vr) d V 

Generalized Boltzmann Equation 

With the aid of these means the discretized version of the generalized Boltzmann 
equation (jBl|) reads 



144W 2 

m,l 



-E'kn -Spm -E qi -Ekn - Pm - qi 

x ^2A cff (t, k n - q b E kn - £ qi ) + Acfr (t, k n - p m , £ kn + E Pm 

X <$(-Ek n + E pm - E ni - ^k„-p m -qij 

x {n kn n Pm (1 + n qi ) (1 + nk n - Pm - qi ) 

- (! + ™k n ) (1 + ™ P J ^q,^k n -p m - qi 
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The discretized 5 function has become a Kronecker 5, i.e. it is 1, if its argument 
vanishes and otherwise. The discretized energy is given by 



We see that due to its character as an integro-differential equation the dis- 
cretized Boltzmann equation becomes a system of coupled ordinary differen- 
tial equations. We recall that the generalized Boltzmann equation stems from the 
kinetic equations for the symmetric propagator and the spectral function which 
were supposed to be valid in the intermediate-time drifting regime, where the evo- 
lution of the system is sufficiently smooth. Of course, this property is inherited 
by the generalized Boltzmann equation, which means that the particle number 
density changes only slightly on its approach to the Bose-Einstein distribution 
function. Therefore we propose to use a fourth order Runge-Kutta method with 
adaptive stepsize control to solve the Boltzmann equation numerically [2H] • Due 
to its adaptive stepsize control this method exploits the smooth behaviour of the 
particle number density and will increase the stepsize very rapidly which allows 
for an efficient approach to the equilibrium. 

The gap equation (|3.4.1j) for the effective mass has to be solved after each 
time step self-consistently. This can be done with a Newton-Raphson method 
38 . Unfortunately this method is not applicable to solve Eqs. (jB3|) to (jB6|) for 
the coupling. Here, one is restricted to an iterative approach to A c ff. 

Examining the collision sum on the right hand side of the discretized Boltz- 
mann equation more closely, we find that the complexity to compute this sum 
explicitely is of order 0(N®). Accordingly, this would make the numerical so- 
lution of the generalized Boltzmann equation an order O (N%) problem, which 
in practice makes the approach to the infinite volume limit impossible. To re- 
duce the complexity of the computation of the collision sum one therefore has 
to use stochastic (Monte-Carlo) methods to estimate its value. Such stochastic 
methods were originally invented to compute higher dimensional integrals by ex- 
ploiting the law of large numbers, namely that for a function / defined in a higher 
dimensional volume V and R random variables X\, . . . ,xr £ V, we have 




and the effective mass satisfies 
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Therefore we can estimate the integral by 

where the error in the estimation is of order 0(l/\/~R). It is straightforward to 
apply this estimation to the collision sum. For randomly chosen m 1; m 2 , m 3 , / l5 1 2 
one performs the sum over / 3 to eliminate the energy Kronecker 5. As we still 
have Ng equations but only one sum left, we have reduced the complexity for 
solving the generalized Boltzmann equation to order O(Nf). Nevertheless, we 
would like to mention that in giving only the order of the complexity we have 
neglected the constants. Especially R may become very large (between 10 3 and 
10 5 ) such that for small volumes there is only little increase in the efficiency of 
our numerics. For example, for N s = 32 and R = 10 4 one saves a factor of about 
3 ■ 10 3 for every collision sum that has to be computed. 

In order to solve the Boltzmann equation, one starts with an initial particle 
number density, i.e. for each k n one specifies a certain value n(0,k n ). For this 
given density, one then has to compute n (ti, k n ) for every momentum mode k n 
after the first timestep, thus obtaining the particle number density for this time. 
This procedure has to be repeated until the particle number density does not 
change anymore up to some given accuracy, such that one can assume that the 
system has reached its equilibrium state. The temperature of the system can 
then be obtained by fitting the particle number density to the Bose-Einstein 
distribution function. 

Kinetic Equations 

It is our aim to trace the evolution of the propagator for some given far-from- 
equilibrium initial conditions with the aid of the full evolution equations (|K1|) 
to (|K6|) until we find that every momentum mode has entered the intermediate- 
time drifting regime. Having arrived at this point, we would like to switch to 
the kinetic equations in order to observe thermalization in a more efficient way. 
The first step, namely the numerical implementation to solve the full evolution 
equations for a spatially homogeneous system, is described in detail in Ref. 0. 
The remainder of this section is devoted to the description of the second and the 
third step, namely the transition from the full evolution equations to the kinetic 
equations as well as the numerical implementation to solve the latter ones. These 
two steps also are described for a spatially homogeneous system. 

Consider Fig. [T] in the introduction and suppose that for some time t all mo- 
mentum modes have entered the intermediate-time drifting regime. For example 
for t = 100, we see that this is the case for at least all momentum modes shown 
in Fig. Then we could choose to to be the time for the transition to the kinetic 
equations. Of course, the larger to the more accurately the kinetic equations 
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describe the subsequent approach to thermal equilibrium. The transition itself is 
done in the following way: Taking the effective memory loss shown in Fig. |2 into 
account, one can introduce a cut-off at the lower limits of the memory integrals: 




X 0-S 



S denotes the size of the storage where we save the values of the propagator for 
previous times. Let a t be the stepsize in time direction and N t the number of 
values stored for each quantity at previous times, then 

S = N t a t . 

In order to perform the transition at the time to we have to compute the evolution 
of the propagator with the full evolution equations up to times x° = y° = to + 
(S/2), where we have stored the values of the propagator for all times 

(x°, y°) e jto - |, t - | + a t , . . . , t Q + | 

For fixed center time to we see that this set contains N t + 1 pairs (x°,y°) with 
different values for s°. Hence, the energy dependence of the Wigner transformed 
propagator is also encoded in N t + 1 different values of the energy and we can 
perform the Wigner transformation according to 

F(x°,y°,k n ) — F(to,s° m ,k n ) 

N t 

— ► F (t , k n , io m ) = a t ^ exp (^ m s°) F (t , k n , s°) , 

j=0 

where m G {0, . . . , N t } and 

2tt 

u m = -g-rn . 

This transformation has to be performed for every spatial momentum mode k n 
and every of the functions F, g, £p, S e , Ip, I s , Gr, S r and Ir. Especially for 
the retarded quantities it is crucial to use the above procedure in order to avoid 
the convolutions with the Fourier transformed 9 function. 

After having performed the transition from the full evolution equations to the 
kinetic equations, one can then use the standard lattice discretization described 
at the beginning of this section to discretize the kinetic equations and solve them 
numerically using standard techniques. Exactly as for the generalized Boltzmann 
equation, we suggest to use a fourth order Runge-Kutta method with adaptive 
stepsize control to advance in time. Exploiting the smooth behaviour of the 
subsequent evolution of the propagator this method will increase the stepsize 
rapidly and thus allow for an efficient approach to thermal equilibrium. 
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Starting from the 1 /iV-expansion of the 2PI effective action to next-to-leading 
order, we derived kinetic equations for a systematic nonperturbative description 
of the evolution of 0(A r )-symmetric real scalar quantum fields out of thermal 
equilibrium. 

The techniques presented here have several generic advantages. First of all, 
evolution equations derived from the effective action are certain to be time- 
reversal invariant, and the corresponding kinetic equations are known to possess 
exact rather than approximate conservation laws on the level of expectation val- 
ues. Additionally, the 1/iV-expansion provides a nonperturbative description of 
quantum fields out of equilibrium. As a consequence the range of applicability of 
our kinetic equations is not limited by the question of the validity of a loop expan- 
sion, but instead is restricted only by the gradient expansion employed in their 
derivation. In situations where this gradient expansion is justified, our kinetic 
equations describe the evolution of quantum fields out of equilibrium without fur- 
ther approximations. In particular we showed that our kinetic equations remain 
valid in the case of nonperturbatively large fluctuations. 

Furthermore, we showed how one can get from our kinetic equations to a gen- 
eralized Boltzmann equation for quasi-particles. While the three-loop approx- 
imation of the 2PI effective action comprises the classical Boltzmann equation 
fHllS], the 1/iV-expansion of the 2PI effective action to next-to-leading order 
leads to a generalized Boltzmann equation which includes an effective mass and 
an effective coupling. 

Finally, this diploma thesis can serve as a starting point for further research. 
In very much the same way as we obtained our kinetic equations in this work, one 
should be able to derive kinetic equations in the case of broken symmetries where 
the field expectation value does not vanish. This would be important to describe 
the phenomenon of parametric resonance for large times. Apart from that, as 
already indicated in the introduction, our kinetic equations are expected to be 
valid near a second-order phase transition. Additionally, we mentioned several 
times that we expect our generalized Boltzmann equation to be valid even for 
large particle number densities, where one cannot neglect collisions between three 
or more particles. Of course, these expectations have to be tested. Further aims 
are the application of these methods to fermionic and gauge theories. 
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